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ABSTRACT 

Pixel microlensing, i.e. gravitational microlensing of unresolved stars, can be used to explore distant 
stellar systems, and as a bonus may be able to detect extragalactic planets. In these studies, binary-lens 
events with multiple high-magnification peaks are crucial. Considering only those events which exhibit 
caustic crossings, we estimate the fraction of binary events in several example pixel microlensing surveys 
and compare them to the fraction of binary events in a classical survey with resolved stars. We find a 
considerable enhancement of the relative rate of binary events in pixel microlensing surveys, relative to 
surveys with resolved sources. We calculate the rate distribution of binary events with respect to the 
time between caustic crossings. We consider possible surveys of M31 with ground-based telescopes and 
of M87 with HST and NGST. For the latter, a pixel microlensing survey taking one image a day may 
observe of order one dozen binary events per month. 

Subject headings: binaries:general — gravitational lensing — planetary systems 



1. INTRODUCTION 

Gravitational microlensing has been a powerful probe 
of massive dark objects in the galactic halo (Paczyhski 
1986; Griest 1991; Alcock et al. 1995). If the lens is 
a binary, a pair of caustic crossings may be observed. 
The caustics of the lens system are curves on which the 
magnification is formally infinite, and quite large in prac- 
tice. In addition to caustic crossing events, binary lenses 
can give smaller perturbations to a standard microlensing 
lightcurve. Such perturbations are in principle detectable 
(Mao & di Stefano 1995; di Stefano & Perna 1997). Bi- 
nary events yield more information about the source-lens 
system than single-lens events, and it is even possible to 
image the surface of the source star due to the large mag- 
nification. Caustic crossing events have been seen towards 
the galactic bulge and the Small and Large Magellanic 
Clouds (Alcock et al. 1997; 2000a; Afonso et al. 1998; 
Albrow et al. 1995; 1999; Udalski et al. 2000; Afonso et 
al. 2000). Mao & Paczyhski (1991) studied caustic cross- 
ing events in classical microlensing. They came to the 
conclusion that approximately 7% of microlensing events 
towards the galactic bulge would exhibit caustic crossings, 
with the assumption of perfect time sampling. 

The pixel technique extends the feasibility of microlens- 
ing experiments by removing the requirement that indi- 
vidual stars be resolved (Grotts 1992, Baillon et al. 1993). 
The possibihty of detecting caustic crossings due to bi- 
nary lenses in pixel events is intriguing, as it allows us 
to study stellar populations at great distances. Several 
groups have successfully used the pixel, or difference im- 
age analysis technique to detect candidate microlensing 
events towards M31 (Grotts & Tomaney 1996; Ansari et 
al. 1999). The MACHO collaboration has also used this 
technique to identify candidates towards the LMC, in ad- 
dition to the classical events (Alcock et al. 1999). This 



technique has a great future as more telescope resources 
are brought to bear, both on the ground and in space. 
Space based surveys will allow the reach of microlensing 
surveys to extend at least as far as the Virgo cluster, at a 
distance of about 16 Mpc. 

In this paper, we extend the Mao & Paczyhski (1991) 
analysis to pixel microlensing surveys of stellar systems. 
In pixel microlensing, detected single-lens events tend to 
have high magnifications, so their rate is suppressed by 
a factor of 1/A, where A is the magnification required 
for detection. However, caustic crossing events always ex- 
hibit large magnifications at the caustic crossings. This 
enhances the rate of caustic crossing events relative to 
single-lens events in pixel microlensing. On the other 
hand, finite source effects will limit the maximum mag- 
nification during a caustic crossing, rendering some clas- 
sically detectable caustic crossing events undetectable in 
pixel microlensing. In this paper we show that the larger 
region of high magnification is the stronger effect. As a re- 
sult, caustic crossing events in pixel microlensing represent 
fractions of 10-15% of microlensing events, compared with 
the 2-3% fraction we find for observable caustic crossing 
events with resolved sources. Finer time sampling would 
increase both fractions, and the latter would then agree 
with the Mao & Paczyhski (1991) result. 

Our primary goal with this work is to show that 
microlensing events exhibiting caustic crossings are de- 
tectable in pixel microlensing surveys. To this end, we 
estimate the rate of caustic crossing events given some 
basic assumptions. Less important to our purposes are 
considerations of accurately determining the properties of 
the source-lens system, which we leave to future work. 
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2. THE TWO POINT MASS GRAVITATIONAL LENS 

We begin with a theoretical discussion of the two point 
mass gravitational lens, the appropriate model of a bi- 
nary star system. This model has been studied exten- 
sively (Schneider & WeiB 1986; Erdl & Schneider 1993). 
We quote some relevant results. We take the masses of 
the lenses to be Mi and M2, with M = Mi + M2, and 
we define dimensionless masses ^1^2 — such that 

Ml + M2 = 1- Furthermore, we define q = M2/M1, with 
< g < 1 without loss of generality. We take the stan- 
dard notation of distances to the source Dg = L, deflector 
Dd = xL, and deflector-source distance D^s = (1 ~ x)L. 
We define the Einstein radius and the Einstein angle for 
the system as (Einstein 1936) 



Re = \l ^^^^x(l 



x), 



Lx 



(1) 



We adopt a complex parameterization (Witt 1990) of 
the lens system. We introduce complex angular coordi- 
nates on the plane of the sky, z = {O^ + i9y)/0E- Given 
two lenses at angular positions zi and Z2, a source at C, 
will have images at the solutions z of the lens equation 
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Here a is the defiection angle, and the bars indicate com- 
plex conjugation. We will call the plane described by C, 
the source plane, and the plane described by z the im- 
age plane. The variables C and z are angles expressed 
in units of the Einstein angle 9e- A source at angular 
position C is at a transverse distance L9e\C.\ — Re\C,\/x 
from the line of sight through the origin of the coordinate 
system. Likewise, a lens at angular position z is at a dis- 
tance xL0e\z\ = Re\z\ from the line of sight through the 
origin. We are interested in the apparent physical size of 
the caustic curves, as this determines the allowed physical 
trajectories of the source which produce caustic crossings. 
The caustic curve, when considered to be in the physical 
source plane, is given by LOeC — ReC/^- 

In this formalism, the Jacobian of C{z) in eq. (0) is 
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dz 



(3) 



The magnification of an image at z is given by l/|J(z)|, 
with the sign of J{z) giving the parity of the image. The 
total magnification is the sum of the individual magnifica- 
tions of the images. 

The points z where the Jacobian vanishes, J{z) ~ 0, 
trace out the critical curves. The inverse images C,i.z) of 
the critical curves are the caustic curves. When the source 
is on a caustic, and so its images are on the critical curve, 
the magnification is formally infinite. 

Let d be the angular distance between the two lenses 
in units of the Einstein angle Oe- According to the value 
of d, there are three regimes in a binary lens system: the 
so-called close, intermediate, and wide binaries. A close 
binary exhibits three caustics, one with four cusps on the 
line of the lenses and two with three cusps off the line. The 
intermediate binary exhibits one six-cusped caustic, and 
the wide binary exhibits two four-cusped caustics on the 



lens line. The separation points c?ci and diw between close 
and intermediate and between intermediate and wide bina- 
ries are as follows (Erdl & Schneider 1993; Rhie & Bennett 
1999), 
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In the equal mass case, fJ^i = fJ-2 = 1/2, we have c?iw = 2 
and dci = l/\/2- As /ii/12 ^ 0, the intermediate binary 
vanishes, and the transition between close and wide is at 
dci = diw = 1- 

The critical curves are simply parameterized by (Witt 
1990) 



At each 0, there are four roots of the resulting quartic 
equation in z. Each root 2^(0) (i = 1,...,4) describes a 
portion of the critical curve. Extending the parameter 
range to < < Stt, these portions can be patched to- 
gether so that they follow each other continuously for every 
disjoint piece of the caustic. In this way, we achieve one 
"continuous" parameterization of the caustic, z = z{(f)). 
The continuous parameterization is explicitly given by 
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Here, fi gives the degree of asymmetry of the lens, b is 
something like a dipole moment about the origin, and the 
other parameters have no obvious simple interpretation. 
The parameterization of the caustics, C = C(0)j follows 
trivially from that of the critical curves, z — z{(f)), using 
the lens equation, eq. (||). 

Lastly, we expand the parameterizations of the critical 
curves and caustics in the limits where d — *■ 0, c». Defining 
u = exp(?(/)/2) with < < An, we find as d ^ dci: 
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with 0((i ^) errors. 



3. RATE OF CAUSTIC CROSSING 

In computing the rate of single-lens events, the Ein- 
stein angle is used as a "cross section" for microlensing. 
Any source which passes inside the Einstein ring is said to 
be microlensed. Griest (1991) has shown that the rate of 
single-lens events depends linearly on the threshold mini- 
mum angular separation between source and lens, imply- 
ing that the dependence is on the angular diameter rather 
than on the solid angle, which would imply a quadratic 
dependence on the impact parameter. In contrast, the op- 
tical depth is a measure of the probability that a given 
star is microlensed at a given instant. This is proportional 
to the area of the Einstein ring. 

To compute the rate of caustic crossing events, we com- 
pute the "angular diameter" of the caustic structure in 
direct analogy to the single-lens case. We will call it the 
"angular width" of the caustic. We will appropriately av- 
erage this "angular width" over the distribution of binary 
systems. The optical depth to caustic crossing would simi- 
larly depend on the solid angle enclosed by the caustic, but 
the rate of caustic crossing events depends on the angular 
width, as we now show. 

In a caustic crossing event, the source crosses the caus- 
tic an even number of times, half entering and half exiting 
the interior of the caustic. Consider the rate Eg at which 
the sources enter the region limited by the caustic. The 
number of caustic crossings per unit time, counting each 
caustic peak in the lightcurve, is then 2Eo. If the orien- 
tation of the caustic is random, as is the case for binary 
systems whose separation has no preferred direction in the 
sky, we find 

To^<S>w, (19) 

where $ is the angle-averaged number flux of the sources, 
and w is the "width of the caustic," deflned in terms of 
the caustic length £ by 

w = -. (20) 

TT 

We take w and £ in units of and $ in the corresponding 
units. 

Eqs. (|l|) and are proven as follows. The flux of 
sources whose proper motion is in a direction forming an 
angle ip with the binary separation is ^dip/{27T). Consider 
an infinitesimal piece of the caustic d£, whose interior nor- 
mal makes an angle 9 with the fixed axis. The number of 
sources entering the caustic through d£ per unit time is 



di; 
2^ 



d£ 

cos{^P - e)d£ = ^— ^<i>dw. (21) 



The limits of integration restrict the flow of sources to 
those coming from one side of the caustic only. Integra- 
tion of eq. ( pi| ) along the caustic gives eqs. ( |l9| ) and (20). 



Eqs. ( |l9| ) and (20) apply to any closed curve in the plane. 
For example, consider the rate of single lens events. The 
critical curve of a single lens is a circle, namely the Ein- 
stein ring. Its circumference is its length £, and eq. (20) 
implies that w is the diameter of the circle. Then eq. (19) 
states the obvious fact that the number of sources entering 
the circle equals the flux times the circle diameter. 

This method of computing the cross section of a closed 
curve counts all entries to the curve. A curve that is not 
convex, such as a caustic curve, can have multiple entries. 
The rate Tq counts all of these entries. For example, a 
binary lens event with 4 caustic crossings is counted twice 
in Tq. In other words, the rate Fq is half the number of 
caustic peaks per unit time. 

4. BINARY SYSTEMS 

We now discuss the observed physical parameters of bi- 
nary systems relevant to our calculation. It is well known 
that a large fraction of stars have companions, perhaps 
the majority. A well-known rule of thumb is that the dis- 
tribution of periods P for binary systems is constant in 
InP, with roughly 10% of binaries in each decade from 
one-third day to ten million years. Less well known is 
the correlation between the masses of the bound stars. At 
high masses, a correlation is seen, but at low masses, the 
data are inconclusive. 

We proceed to express the separation of the binary pair 
in terms of its period and total mass. With semimajor 
axis o, Kepler's law reads 



f— V 
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(22) 



The Einstein radius is 



Re 
AU 



= B 




(23) 



with B = 0.28538\/a;(l - x)10^/i°, and D is the distance 
modulus to the source. Assuming that the projected sep- 
aration is at maximum, the most probable situation, we 
find the relation 



2/3 



M 



-1/6 



(24) 



In the range P/yr = [10^^, 10^], the distribution in P is 
simply dN/dlogP = 1/10, and is correctly normalized. 
Fixing all other quantities, appropriate for the rate calcu- 
lation we will do, we thus find the distribution of binaries 
in d, dN/dlogd = 3/20. 

There is disagreement about the distribution in q of 
binary systems, ranging from a linear rise cx 2.6 + 2.9q 
(Mazeh et al. 1992) to a power law dechne cx q~^-^ (Pa- 
tience et al. 1998), with other measurements falling be- 
tween (Duquennoy & Mayor 1991; Trimble 1990). In the 
next section we will find that the effects of the q distribu- 
tion of binaries are mild. 
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5. SIZE OF CAUSTIC STRUCTURE 

We evaluate the mean width of the caustics w for the 

binary systems described in the previous section. In the 
parametric representation of the caustic C, = C(0)) its 
length in units of the Einstein angle is 



d<l) 



dcj). (25) 



We have computed the length of the caustic structure 
for various values of q. The caustic length is not as sen- 
sitive to the q distribution as it is to the lens separation 
d, as we show in figure 1. There we plot the width of the 
caustic as a function of d, averaging over four distribu- 
tions in q: a linear rise oc 2.6 + 2.9g, flat, flat in \nq, and 
a power law q~^'^ (in the last two distributions we take 
q > 0.1). The difference is only a few per cent. The peaks 
occur at the average values of the close; intermediate and 
intermediate-wide binary separations dci and diw- 




lens separation d 

Fig 1. — Average width w of caustic structure as a function 

of lens separation d, both in units of the Einstein angle 
0E- We have averaged over several distributions in the bi- 
nary mass ratio q, and also included the case q = 10~^, of 
interest for planetary binaries. We have also plotted the 
mean chord length, averaged over the flat q distribution. 

For small and large separations d, we have found the 
leading behavior of the mean caustic width. In both cases 
the important caustics are approximately square. For 
d < dci- the linear dimension of the two additional tri- 
angular caustics is suppressed relative to the single, cen- 
tral square caustic by a factor of d. The mean width 
is w = 3^2(1 



/i )/7r in the case d <§; dci, and it is 



w = 6(1 — |I)^/2(T+Jl)/{'Kd'^) when d > diw Thus, we 
see that the caustic structure is largest when d is of order 
unity. Averaging over the four q distributions mentioned 
above (namely 2.6 -|- 2.9g, flat, 1/q, and q~^'^, in this or- 
der) we find w/d'^ = 0.78, 0.74, 0.68, and 0.64 in the case 
d <C dci, and wd"^ = 2.0, 1.85, 1.63, and 1.51 in the case 
d > diw 

Now that we have computed the cross section for caustic 
crossing events, we can compare with that of single-lens 



events on resolved stars, whose cross section in Einstein 
units is 2. The cross section averaged over angles is £/tt, 
which needs to be averaged over d and q. We find the ra- 
tio of caustic crossing to single-lens events to be approx- 
imately 6.3%, multiplied by the fraction of star systems 
that are binaries. This is in quite good agreement with 
the results of Mao and Paczyhski (1991). 

In addition to the mean width of the caustic structure, 
important for determining the cross section for caustic 
crossing events, we would like to know the mean length 
of chords through the caustic structure. This quantity is 
directly related to the mean time between caustic cross- 
ings. We proceed as follows. Let w{9) be the projected 
width of the caustic structure at a fixed orientation 0, and 
let C(9, b) be the chord length at orientation 9 and impact 
parameter b. The area of the caustic structure is clearly 



dbC{e,b) = w{e)C{e), 



(26) 



where now C(9) is the mean chord length at orientation 
9. The area is of course independent of 9. We want the 

mean chord length, weighted by the width of the caustic. 
The mean chord length over all orientations is now just 



C = A 



d9 1 w(9) 



A - ^' (27) 
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the area divided by the mean cross section. Note that 
on the occasions where there are multiple pairs of caustic 
crossings, the chord is counted as multiple separate pieces. 

We will need the area and projected diameter of the 
caustic structure. In vector notation, we have dA = 
{fx df)/2, which is easily transformed to the complex 
parameterization, 



A = 



/■Stt 



Im C 



(28) 



j\m (e-^-^/^c) Re (^e-*^/2 



dz 



The projected diameter along the real axis, at orientation 
9, is simply 



w{9) 



1 



Sir 



|dRe(e'^C)| 



(29) 



cos I 6* + 1 Rc 



-i<p/2 



dz 



Note that if this expression is averaged over orientations 
9, we recover the result for the mean width of the caustic 
structure being proportional to the length of the caustic. 

We have now assembled all of the necessary pieces for 
calculating the mean length of chords through the caustic 
curves. As for the mean diameter, we compute this quan- 
tity as a function of the lens spacing d, averaging over a 
distribution in q, which we take to be constant in q. The 
mean chord length is plotted in figure 1. 
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We arc furthermore interested in the fuU distribution of 
chord lengths for a given lens configuration. As is done 
by Han, Park & Lee (2000), we have performed a Monte 
Carlo simulation to find this distribution. The results are 
shown in figure 2, with the exact mean shown as a short 
dashed line and the median shown as a long dashed line. 
From this simulation, it seems that using the mean chord 
length may not be an adequate prescription for designing a 
microlensing sm vey, since the distribution of chord lengths 
is quite complicated. However, wc will find that the mean 
chord length prescription is usually quite good. 
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Fig 2. — Distributions of chord lengths for several choices 
of d and q. The exact means of these distributions are 
indicated by the short dashed lines, while the medians are 
shown by the long dashed lines. The small panels show the 
shapes of the caustic curves, with the masses at x = ±d/2 
and the larger mass on the right. The y-axis is in arbitrary 
units. 

6. LIGHTCURVES AND FINITE SIZE EFFECTS 

We illustrate two lightcurves for caustic crossing events 
in figure 3. Other examples can be found in the litera- 
ture cited above. In the upper panels wc show the caustic 
structure together with the source trajectory in the plane 
{9x,0y). We place the lenses along the ^^-axis at ±d/2. 
In the lower panels wc plot the lightcurves as functions 
of 0x/9e, which for a uniformly moving source is related 
linearly to time. 

In the interior neighborhood of a caustic, the magnifi- 
cation diverges as K/ where y is the angular distance 
perpendicular to the caustic, in units of the Einstein angle 
(Chang & Refsdal 1979, Kayser & Witt 1989). The flux 
factor K depends on the lens map through the length of 
the tangent vector to the caustic, T = \T^\, as K = y/2/T. 
In the complex parameterization, the tangent vectors to 
the critical curve and caustic are (Witt 1990) 



-2i 



.da d'^a 
dz dz'^ ' 



oz 



(30) 



Defining the Einstein time as usual ts = 2Re/v, and 
taking the angle to the normal as 6, the lightcurve near 



the caustic crossing at t = to is, in the interior side of the 
caustic, 



tE 



FqK 
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(31) 



with the + sign (— sign) when entering (exiting) the caus- 
tic. In pixel lensing Fq is unknown (but see the next para- 
graph for a possibility to determine it). However, there 
are two caustic crossings, and if the timescales of the two 
are compared, the ratio of the geometric factors K/y/cos6 
at each of the two caustic crossings can be measured. This 
may provide some insight into the exact parameters of in- 
dividual events. 
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Fig 3. Example source trajectories and lightcurves for 
caustic crossing events. The lenses are marked with 
crosses, the more massive on the right. 

We may gain useful information by measuring the du- 
ration of the caustic crossing. This is the time over which 
the disk of the source star intersects the caustic, and is 
given simply by 



tE 



2Lx9, 



9e cos 6 V cos 5 



(32) 



For star-star lensing in a distant galaxy, x k, 1, L is known, 

and V and cos 5 are taken from known distributions. This 
can give a handle on the angular size of the source star. 
Furthermore, if the caustic crossing is seen in multiple 
wavebands, the color(s) of the source star can be deter- 
mined. This is possible because the magnification is the 
same in all wavebands, thus comparing the flux increase 
in two bands is the same as comparing the actual (unmea- 
sured) source flux in those bands. With the color and the 
angular size, a rough estimate of the absolute flux Fq can 
be derived. It may also be possible to image the surface of 
these extragalactic stars as the caustic passes over them, 
as has been done towards the SMC (Afonso et al. 1998; 
Albrow et al. 1999). However, this would require intense 
monitoring to cover the actual caustic crossing. It is not 
the purpose of this paper to discuss this possibility further. 
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The magnification during a caustic crossing is critically 
important to a calculation of the event rate. While the 
magnification of a point-like source is infinite, that of a 
real source is finite. Since only magnifications higher than 
a threshold magnification can be observed, it is crucial to 
determine the fraction f{A>Ao) of caustic crossing events 
that exceed a fixed magnification Aq. The rate of caustic 
crossings entering the caustic with magnification A > Aq 
is then 

T^fiA>Ao)To. (33) 

We now determine f{A>Ao). 

The maximum magnification of a finite source of angular 
radius 6s is given by (Schneider & Weii3 1987) 



A„ 



(34) 



where fs is a form factor depending on the luminosity 
profile of the source. For a uniform disk, fs — 1.39, while 
a limb-darkened profile ~ a/I — r^/i?^ gives fs — 1.47 
(Kayser & Witt 1989), thus we see that the profile is rel- 
atively unimportant for our purposes. 

AmnK in eq. (p3) is the highest magnification when the 
source crosses tEe caustic. It follows that the fraction 
of events that cross a given caustic with magnification 
A > Aq is equal to the fraction of length of the caustic 
in which K>Ko = Aov/^ry^//,, 



fiA>Ao) 



w{K>Ko) 



(35) 



Taking a distribution of mass ratios constant in q, we plot 
w{K>Ka) as a function of d for several values of Kq in 
figure 4. 
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Fig 4. — Effective width of the caustic structure w{K>Kq). 
We choose values equally spaced in logKa. The curves de- 
crease in magnitude as increases. 

7. EVENT RATE 

We now compute the rate distribution of caustic cross- 
ing events. We let the velocity distribution of the source 



and lens population be Maxwellian with characteristic ve- 
locity Vc- We denote the transverse velocity of the lens 
vt , and the transverse velocity of the line-of-sight vi . We 
denote the total mass of the lens M, and the apparent 
magnitude of the source m. We define the mass function 
£,{M) ~ dN/dM normalized to unity, and the luminosity 
function (j){m) normalized to the surface brightness. We 
denote the distribution in d as n(c?) = 3/(20(iln 10). The 
mean width of the caustic structure in Einstein units is w, 
and we choose a distribution constant in q. 

We will compute the rate of caustic crossing events with 
respect to the time between caustic crossings tc in two 
ways. First we will use the mean length of chords to de- 
termine this timescale, which is not correct event by event, 
but may be acceptable for a rate distribution. If we define 
the effective time between caustic crossings 



Re{C) 



Vt 



(36) 



where (C) = {C{d)) is the length of the caustic chord 
averaged over the q distribution, we can write the rate 
distribution as 



dT 



dt 



2Lvl I dm(f){m) / dM 



e(M) 
M 



dxp'{x) (37) 



ddn(d).^e-(--)^/„(2.,)^4^. 

with dniin given below, lo — RE(x){C{d)) j {vct(c)) and 
r/ = vi/vc- Also, Iq(x) = /o(a;)e~^, where Iq is a mod- 
ified Bessel func tion of the second kind. We note that 
/o(x > 1) w l/%/2^ and io(0) = 1. 

Secondly, we will use the actual caustic crossing 
timescale, at the cost of doing an exhaustive Monte Carlo 
of chord lengths. Here we define tc as in eq. (p^), except 
that now C is the actual chord length. The rate distribu- 
tion is then 

^^2Lvj [ dmcl){m) [ dM [ dxp'{x) (38) 



dtr 



M 



I dt;w3e-(^-'')'/o(2r;i;) Weff f 



Vctc 

Re 



V, K > Kc 



with V = Vt/vc, WcS being the effective width as a function 
of chord length. 



Weft (C, Ko 



dN , 



(39) 



e,q,d 



w is the projected diameter given in eq. (^9|), and dN/dC 
is the probability distribution of chord lengths at a given 
9, q, d, computed by Monte Carlo and normalized to unity 
at Kq — 0. We illustrate the function for several val- 
ues of Kq in figure 5. 

To fix the minimum flux factor Kq, we adopt the con- 
vention that the signal to noise totaled throughout the 
approximate lightcurve eq. (|3l|), be greater than Qq = 7. 
We assume that samples are taken daily, (except for the 
reference case of the Milky Way Bulge, monitored con- 
tinuously). The lightcurve is truncated first by assuming 
the first sample is taken at one half day away from the 
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caustic crosHiiig, and then by requiring that the steUar di- 
ameter is less than the; length of the chord. We derive 
the signal to noise as follows. Assume that the telescope 
collects photons from a source of magnitude m at a rate 
5*0 10"° '''". An integration of tint per observation is taken. 
The background galaxy light falling on the pixel has a 
magnitude fi. This gives a signal to noise per sample of 
Q = v%t~;iO-°-4('"-^/2)(^ _ .^l^ere we have only 

considered photon noise. Ground based microlensing sur- 
veys have achieved noise levels within a factor of two of 
the photon noise, and we expect space based surveys to 
approach the photon noise limit. Thus there is a simple 
relation between the minimum allowed Q = Qo and the 
minimum flux factor Kq. 



on y and L and not on s), we find 




0.5 1 
chord length C = 2t^/t^ 

Fig 5. — Effective width for caustic crossing events. Ten 
million trajectories through caustic structures have been 
generated for various values of d, q, and 9, and the his- 
tograms appropriately weighted for several values of Kq. 
The Weff histograms have been multiplied by C for clarity. 

The minimum d allowed is determined by fixing 
the period at its minimum, 10^'^ yr, giving dmin = 
O.O2(M/M0)-i/6B(a;)-i. In practice, fixing rf^in = 0.1 
is fine, since d lower than this gives a very small cross sec- 
tion, and for the mass we are concerned with, it gives a 
minimum lens-source distance of 10 pc, which is negligi- 
ble. In principle a maximum d can be determined by the 
binding of the binary, but it is never important observa- 
tionally. 

The density p'{x) is the effective density of source-lens 
systems at separation x, and can be determined by inte- 
grating the source and lens positions along the line of sight 
through the target system. Assuming a radial density pro- 
file of stars p{r) truncated at a radius R, we calculate the 
column density along a line of sight with impact parameter 
b, 

rVR^ 



N 



-I 



ds. 



(40) 



The density along the line of sight is used as a probability 
distribution for sources. Using y = 1 — x, and assuming 
L ^ s (which furthermore implies that Re only depends 



p'{l-y)= / -p(V&2 + s2) X (41) 



p (^62 + 8"^ + 2Lys + L'^y'^^ ds. 

Finally, we comment on the rate of binary events with 
resolved sources. In this case there is no maximum x, as 
any caustic crossing event can be detected. Thus the in- 
tegral over the luminosity function (/)(m) is trivial. The 
remaining triple integral is simply multiplied by the num- 
ber of monitored stars. 

8. BACKGROUNDS AND BLENDING 

In any microlensing survey, the issue of backgrounds 
is crucial. In considering binary microlensing events, the 
lightcurves are basically characterized by a double spike. 
Such events in isolation should be exceedingly rare. The 
most pernicious of backgrounds for single lens events is the 
mira variables (Crotts et al. 2000), which have a very long 
period, but as they are characterized by a single spike, 
we are unconcerned with them. However, if one of the 
two caustic crossings is missed, the characterization of the 
event as microlensing would be very difficult, and only pos- 
sible with a very high signal to noise and/or long follow-up 
period. 
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Fig 6. — Simulated lightcurves in the ACS. The binary 
system is the same as in the left panels of figure 3. We il- 
lustrate lightcurves where the time between caustic cross- 
ings is 3, 5, 10, and 20 days. The brightness is 27.5 mag 
pixel" ^, the source star has an apparent magnitude of 28.5 
in 7, and we have accounted for the ACS PSF. The dot- 
ted curves indicate the theoretical expectation, including 
finite source effects. 

In uncovering variability, we must be concerned with 

the amount of blending of the lensed star's light. We will 
now make a simple estimate of this quantity. We want to 
know first the typical magnitude of the lensed star. Study- 
ing a typical /-band luminosity function, such as that of 
Terndrup, Frogel & Whitford (1990), we see that the slope 
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at the highest fluxes is shaUower than 0.4, meaning that 
the most hght comes from the brightest stars. This means 
that the majority of detectable microlensing events wifl 
be microlensing of very bright stars. In fact, typically the 
rate of such events is expected to peak somewhere between 
the surface brightness fluctuation magnitude (which for 
the Terndrup et al. (1990) luminosity function is roughly 
Mi = —1.5) and the tip of the red giant branch (at 
roughly M/ — —3.5). As a rough estimate for M87, if we 
assume the pixel scale of the Advanced Camera for Sur- 
veys (ACS) on the HST is 0.05", and a surface brightness 
of roughly 21 mag arcsec^^ for most of the field of view, 
the pixel brightness is 27.5 mag pixel^^, while the sur- 
face brightness fluctuation magnitude is roughly 29.5 for 
M87, meaning that each pixel has the light of roughly 6 
SBF magnitude stars. Moreover, we assume a PSF where 
roughly 25% of the light of a star falls on the central pixel. 
This means that if we desire the pixel flux to double, the 
typical required magnification (for Mj = —2.5) is a factor 
of 12. In an 1800 second integration with the ACS at this 
surface brightness, roughly 350 photons are collected, so a 
IOct excess requires only a 50% magnification of the pixel 
flux, or a microlensing magnification that is typically of 
order six. 

In figure 6, we illustrate four representative lightcurves 
of a binary event in M87 observed by the HST ACS, sam- 
pling daily. In all cases the geometry is identical to the 
left panel of figure 3, but with four different timescales be- 
tween caustic crossings of 3, 5, 10, and 20 days. With a 
timescale of three days, though the event is detected at a 
significant signal to noise, it is not convincingly a binary. 
With a timescale of five days or more, it seems more likely 
that binary events can be unambiguous. 

Finite source size effects are not severe for the example 
lightcurves in figure 6. The star of magnitude M/ — —2.5 
has a radius of roughly 65 Rq (Gould 1995), and for a 
solar mass lens system with Dds ~ 5 kpc, the star sub- 
tends roughly one seventh of the angular distance between 
caustics along its path. The dotted curves are the theo- 
retical lightcurves, including finite source effects for a uni- 
form surface brightness star, computed according to Gould 
& Gaucherel (1997). According to eq. (^, the hmiting 
counts should be roughly 475 on the left and 600 on the 
right. As shown, the full calculation is not quite that se- 
vere. That the lightcurves should not be seriously affected 
is also clear from the fact that the angular size of the star 
is small relative to the angular distance between the caus- 
tics. For the two longer events, it becomes possible to see 
the actual caustic crossing happening, though with one 
image a day the accuracy is poor. In principle, one could 
envision an alert system so that more frequent images are 
recorded during a suspected caustic crossing. 



9. EXAMPLE MICROLENSING SURVEYS 

We consider the case of star-star lensing in M87, a giant 
elliptical in the Virgo cluster at a distance of 15.8 Mpc. 
We take the radial profile of stars to be 



with 



1,1 + 0.275 log 



(43) 



(42) 



and po = 3.76 M© pc-^ following Tsai (1993). We adopt 
the /-band luminosity function of Terndrup et al. (1990) 
adjusted to yield a surface brightness fluctuation mag- 
nitude of Mj — —1.5, and a mass function M^{M) oc 
exp(-.4171nM - 0.0886 In^M) with M in solar masses 
(Miller & Scalo 1979). The velocity dispersion is taken to 
be Vc = 360 km s~'^. Typically, stellar radii are Og ~ 10""'^'^, 
while the typical Einstein radii are of order 9e ^ 2x 10^^^ 
for solar mass lenses, a factor of twenty larger. We com- 
pute the rate of caustic crossing events observed in surveys 
using the ACS on the HST, and on a model NGST. We 
compare with the rate of single-lens events (Baltz & Silk 
2000), requiring seven samples detected at 2a, as Crite- 
ria A of Alcock et al. (2000b). We assume that the sensi- 
tivity of the ACS is 4.5 times that of the WFPC2 in the 
I band, meaning that one photon per second is collected 
from a star of magnitude m/ — 25.77, integrated over the 
entire PSF. We model the NGST as having seven times the 
efficiency of WFPC2, but with nine times as much collec- 
tion area. These rates are plotted in the upper panels of 
figure 7. In a monitoring program of daily observations of 
M87 for a period of a month, the Advanced Camera for 
Surveys on HST would detect of order one caustic cross- 
ing event, and the NGST could detect of order one dozen 
caustic crossing events. 

Microlensing in M31, the Andromeda Galaxy, has been 
searched for extensively using the pixel technique (Ansari 
et al. 1997, 1999; Melchior et al. 1998, 1999; Crotts and 
Tomaney 1996; Tomaney and Crotts 1996). We have com- 
puted the rate of caustic crossing events observed in the 
bulge of M31, assuming a CFHT-class telescope and us- 
ing the bulge model of Kent (1989), with a velocity dis- 
persion of Vc = 220 km s^^. The typical star has a ra- 
dius ^ 2 X 10^^^, and the typical Einstein radius is 
~ 4 X 10^^^, again a factor of 20 larger. The fiducial 
star has magnitude m/ — 26.88, giving one photon per 
second integrated over the PSF. Furthermore, we assume 
that the errors are twice the photon counting noise, as 
found by the AGAPE collaboration (Ansari et al. 1997). 
We again compare with the rate of single-lens events, and 
plot the results in figure 7. A CFHT survey of the bulge 
of M31 would detect on the order of one binary event per 
month of observations. 

Finally we compare the rates of single-lens and caus- 
tic crossing events towards the bulge of the Milky Way. 
We take the simple bulge and disk model of Evans (1994), 
with no bar. These rates are plotted in the upper panel of 
figure 7. Note that typically 10^ stars are monitored in the 
galactic bulge. We notice that our rates are comparable to 
the the results of the MACHO project monitoring of the 
galactic bulge (Alcock et al. 1997; 2000a). While we find 
that roughly 6% of events should exhibit caustic crossings 
in this case, the finite stellar radii remove a significant 
fraction of these events as the stellar radius becomes com- 
parable to the distance between the caustics, and we find 
that roughly 2-3% of events will exhibit detectable caustic 
crossings. 

We find that caustic crossing events typically make up 
10%, and as much as 15%, of the total rate of stellar events 
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in pixel microlcnsing surveys, assuming that all stars are 
binaries. This is significantly higher than the relative 
rate in microlensing with resolved sources, typically 2-3%. 
Thus, the larger effective cross section is more important 
than the more severe finite source effects. 

These estimates all assume only self-lensing in the tar- 
get galaxy. Of course there is the possibility that lenses in 
the halo of the Milky Way could contribute. Unless there 
is a significant population of halo lenses, the rate for stel- 
lar lenses in the Milky Way would be substantially smaller 
than the self-lensing rate in the target galaxy, since the 
scale height of the disk is much smaller than the size of 
the spheroids of the target galaxies. 
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Fig 7. — Rate of microlensing events. In all cases the solid 
lines indicate the rate of single lens events as a function 
of the full- width at half maximum timescale ffwhm, the 
dotted lines indicate the rate of caustic crossing events as 
a function of crossing time tc, and the dashed lines in- 
dicate the caustic crossing rate using the mean crossing 
time t(^c} ■ Top panel: M87 events using the Advanced 
Camera for Surveys on HST. Second panel: M87 events 
using the NGST. Third panel: Milky Way bulge events 
for resolved stars. This rate is quoted per monitored star. 
Bottom panel: Events in the bulge of M31. We use the 
improved peak threshold trigger for single lens events in 
M31 and M87. A Miller-Scalo mass function is assumed 
in all systems. 

Prom figure 7 we see that the time between caustic cross- 
ings is typically of order several days. Sampling daily 
would allow a sizeable fraction of the lightcurves to be 
well-measured. We also see that for long timescale events, 
the timescale ti^c) gives a very good approximation to the 
true rate (using tc), but for the shorter events, it is less 
accurate. We notice that the agreement between the two 
is excellent for the NGST and for the M31 bulge. In these 
cases, the required magnification is quite modest, of order 
a factor of two for a star with the surface brightness fluc- 
tuation magnitude. The ACS requires significantly higher 
magnifications, spoiling the excellent agreement between 
the rate with respect to the two timescales tc and i(c)- 
For the Milky Way Bulge, the agreement between the two 



timescales is quite good, though the (quite broad) peaks 
disagree by a factor of about two in timescale. We thus 
find that the timescale ti^c) provides a surprisingly good 
estimate of the rate of events with respect to the time 
between caustic crossings. 

Finally, we comment that with daily sampling, the fine 
details of the caustic crossing will not be well covered. If 
we choose a finer sampling scheme, with the same tele- 
scope resources, the expected number of events will be 
much smaller. We feel that it is important first to identify 
binary lensing events. If this is shown to be possible, then 
more sophisticated sampling strategies should of course be 
considered. 

10. EXTRAGALACTIC PLANETS 

As an aside, we apply the previous analysis to a pixel 
microlensing search of extragalactic planets. 

Discovering planets outside our own solar system has 
been an important scientific goal for many years. The 
planets that have been discovered so far were all found 
by detecting the wobble of the star that they orbit. The 
first detection was of a planet approximately the mass of 
Earth, in orbit around a pulsar (Wolszczan & Frail 1992; 
Wolszczan 1994). More recently, planets have been de- 
tected around main sequence stars (see Marcy, Cochran 
& Mayor 2000 for a review). We would like to extend 
the range of planetary searches to distant galaxies. Gott 
(1981) was the first to point out that gravitational lensing 
might be used to detect low-mass objects in distant galax- 
ies. Here we make another suggestion: pixel microlensing 
surveys may be able to detect Jupiter-mass planets as far 
away as the Virgo cluster, by capitalizing on the unique 
nature of binary microlensing events. 

Previous work has shown that planets might be de- 
tected in microlensing events in the bulge of the Milky 
Way galaxy, or in the Small and Large Magellanic Clouds, 
which are small galaxies in orbit around the Milky Way 
(Mao & Paczynski 1991, Gould A. & Loeb, A. 1992, Gri- 
est, K. & Safizadeh, N. 1998). Stars in the bulge and in 
the Magellanic clouds can be resolved easily, and surveys 
routinely monitor of order ten million stars for microlens- 
ing events. Evidence for a planet orbiting a binary star 
system in the Milky Way bulge has recently been claimed 
in a joint publication of the MPS and GMAN collabora- 
tions (Bennett et al. 1999), although the data can also be 
interpreted as a rotating binary system without a planet 
(Albrow et al. 2000). 

A solar-type planetary system can be described to first 
approximation as a binary object, consisting primarily (in 
mass) of the Sun and the planet Jupiter. In such systems 
q is very small, for example q « 10^''* for the Sun and 
Jupiter. In figure 8, we show two example lightcurves of 
microlensing events, together with the trajectories of the 
source stars relative to the caustic curves. In both cases 
the star has a companion one one-thousandth as massive, 
like the Sun-Jupiter system. 

We have calculated the rate of planetary events observ- 
ing M31 with a telescope like the Canada- France-Hawaii 
telescope (CFHT) on Mauna Kea. We assume that ev- 
ery star has a companion that is one one-thousandth as 
massive, just like the Sun and Jupiter. Furthermore, we 
assume, as is true for known binary stars, that the distri- 
bution of orbital periods is such that ten percent of such 
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systems lie in each decade of period, from a third of a 
day to ten miUion years. This gives a 10% probabiHty 
that a star has a companion between one and five AU, 
in rough agreement with observational findings (Marcy, 
Cochran & Mayor 2000). With a long-term monitoring 
program observing four times daily (using the necessary 
global network of telescopes) , and relaxing the threshold 
to be three sigma in the two images appearing inside the 
caustic, we expect about one caustic crossing planetary 
event every two years. This neglects other types of mi- 
crolensing anomalies, such as those discussed by Covone, 
de Ritis & Marino (1999). 

To increase the chances to detect planetary systems 
in distant galaxies, we require a space telescope such as 
the proposed Next Generation Space Telescope, to be 
launched around 2009. This will be a large (about 7 me- 
ters in diameter) infrared telescope at a Lagrange point of 
the Earth-Sun system, and it will be more than ten times 
as sensitive as the Hubble Space Telescope. 
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Fig 8.— Microlensing events due to a planetary system. We 
illustrate the trajectory of a source star over the caustic 
curves of a lens star and its planet. We then illustrate the 
observed magnification of the source star along its trajec- 
tory. The cross indicates the star's position, whereas the 
planet lies off the plots at —0.5 and 0.6, respectively. The 
plots are in units of the Einstein angle Be, which is the 
characteristic angular scale of a microlensing event. Also 
shown is the star-planet separation d in units of 6e- 



We have calculated the rate of events we might detect 
with the NGST observing the giant elliptical galaxy M87 
in the Virgo cluster, at a distance of 15.8 Mpc. With the 
same assumptions as the M31 calculation, we find that 
an NGST survey of three month's duration, taking four 
images each day, should be able to detect of order four 
planetary systems. We find that such a survey is most 
sensitive to events where the separation between caustic 
crossings is about five days. An alert system for microlens- 
ing events would allow more frequent measurements of the 
light curve during the caustic crossings, with the possibil- 
ity of determining the orbital parameters of the planetary 
system. 

Our results for the rate of caustic crossing planetary 
events indicate that the rates are probably too small for 
this technique to be a feasible method of detecting extra- 
galactic planets. We have made overly generous assump- 
tions, and still the detection probabilities are marginal. 

11. CONCLUSIONS 

We have found that caustic crossing binary events 
should typically make up 10%, and as much as 15%, of 
the total rate of events in pixel microlensing surveys. This 
is significantly higher than the relative rate in microlensing 
with resolved stars, which is typically 2-3%. The enhance- 
ment is due to the larger magnification of caustic crossing 
events with respect to single-lens events. The suppression 
in the count rate expected from finite-source effects is not 
large enough to win over the enhancement from larger re- 
gions of high magnification. 

In examples of pixel lensing surveys, we find that on 
the order of one binary event per month of observation 
should be observable in a CFHT survey of the bulge of 
M31. Also of the order of one caustic crossing event per 
month of daily observations of M87 should be detectable 
with the Advanced Camera for Surveys on HST. The lat- 
ter rate would increase to of order one dozen caustic cross- 
ing events per month with the NGST. Searches for extra- 
galactic planets with the microlensing technique discussed 
in this paper do not seem particularly promising. 
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ful comments and hospitality at Oxford University, where 
this research was begun. E. B. was supported in part by 
grants from NASA and DOE. 
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